Mon. Not. R. Astron. Soc. 000, ??-?? (2001) 



Printed 1 February 2008 



(MN WF&l style file vl.4) 



Cataclysmic variable evolution: AM Her binaries and the 
period gap 



R. F. Webbink 1,2 & D. T. Wickramasinghe 1,2 

1 Department of Astronomy, University of Illinois, Urbana- Champaign 

2 Department of Mathematics, Australian National University, Canberra 



(N 
O 
O 
(N 

u 
< 



Accepted. Received 



in 

m 

o 

O 

Or 

6 



ABSTRACT 

AM Her variables - synchronised magnetic cataclysmic variables (CVs) - exhibit a 
different period distribution from other CVs across the period gap. We show that 
non-AM Her systems may infiltrate the longer-period end of the period gap if they 
are metal-deficient, but that the position and width of the gap in orbital period is 
otherwise insensitive to other binary parameters (excepting the normalisation of the 
braking rate). In AM Her binaries, magnetic braking is reduced as the wind from 
the secondary star may be trapped within the magnetosphere of the white dwarf 
primary. This reduced braking fills the period gap from its short-period end as the 
dipole magnetic moment of the white dwarf increases. The consistency of these models 
with the observed distribution of CVs, both AM Her and non-AM Her type, provides 
compelling evidence supporting magnetic braking as the agent of angular momentum 
loss among long-period CVs, and its disruption as the explanation of the 2 h — 3 h period 
gap among nonmagnetic CVs. 
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1 INTRODUCTION 



The disrupted magnetic braking model has enjoyed consid- 
erable success in explaining the 2 — 3 hour period gap in the 
period distribution of the cataclysmic variables (CVs). The 



in thermal equilibrium throughout its evolution, it would 
evolve continuously through the period gap with no abate- 
ment of mass transfer. 

A subset of the CVs are known to be magnetic and clas 
sified as eith er AM Hers or intermediate polars (IPs 



idea is that the CVs with orbital periods P orb < 10 h terson 1994J ). In the AM Hers the magnetic white dwarf 
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evolve towards shorter periods due to a combination of angu- 
lar momentum losses from magnetic braking due to a mag- 
netised wind originating from the secondary star, and from 
gravitational radiation (Ver bunt & Zwaan (1981); Rappa- 
port, Verbunt & Joss ( 1983| ); see Howell, Nelson & Rappa- 
port (2001) for a recent review). Magnetic braking domi- 
nates at longer periods. As the donor secondary decreases 
in mass, its thermal time-scale increases, until thermal re- 
laxation can no longer keep pace with the mass transfer rate 
being driven by magnetic braking, and it becomes increas- 
ingly bloated with respect to thermal equilibrium. When, 
as the orbital period reaches (~3 h ), the secondary's connec- 
tive envelope penetrates to the stellar center, its large-scale 
magnet ic field loses its radiative anchor and, being buoyant, 
is expelled from the star. This sudden cessation of magnetic 
braking allows the distended secondary to contract more 
rapidly than its Roche lobe, and mass transfer ceases un- 
til angular momentum losses to gravitational radiation have 
reduced the binary separation far enough for the secondary 
to again fill its Roche lobe, at an orbital period of ~2 h . It is 
important to recognize that, were the secondary to remain 
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is in synchronous rotation with the orbit so that its spin 
period P sp m is equal to the orbital period P or b- In the 
IPs, the magnetic white dwarf is asynchronous (typically 
Pspin ~ O.lPorb). Some 65 confirmed AM Hers are known, 
of which 37 have direct ma gnetic field measurements (see 
Wickramasinghe & Ferrario ( |2000| ) for a recent review) . The 
estimated effective dipolar fields range from ~7 — 200 MG. 
The IP class contains fewer confirmed members (26). Indi- 
rect estimates of their field strengths suggest dipolar fields 
in the range ~ 5 — 30 MG. Any theory which purports to 
explain the period distribution of CVs must also account for 
the properties of the magnetic CVs which, as a gro up, com- 
prise some ~20% of all CVs of known orbital period ( Downcs 



et al. 2001) 



The possibility that the magnetic field of the white 
dwarf may affect the orbital evolution of the binary was 
first discussed by Wickramasinghe & Wu ( L994 ) in relation 



to the AM Her type systems. This research was prompted 
by the early discovery of several AM Herculis systems in 
the period gap before the results of the enlarged sample of 
AM Hers from the ROSAT survey became available. Here 
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it was argued that at strong fields, the magnetised wind 
from the secondary may become trapped by the combined 
magnetosphere of the primary and the secondary and result 
in a reduction in the magnetic braking rate. If the braking 
rate was effectively curtailed, systems will evolve towards 
shorter periods almost entirely by gravitational radiation. 
The proposal was that there will be no disruption in mag- 
netic braking, and that no period gap wi l l form. 



Li, Wu & Wickramasinghe (| 1994a| ; [1994b] |1995| ) ex- 
tended the Mestel & Spruit 

dl987| r 

wind theory for single 
stars and developed a theory for evaluating the braking rate 
in synchronised AM Herculis type binary systems. These 
models assume that the secondary star has a dynamo gen- 
erated field with a dipole moment oriented perpendicular 
to the orbital plane, and that the magnetic white dwarf is 
locked in synchronous rotation by magnetostatic torques. 
The dipole moment of the white dwarf is then oriented 
anti-parallel to that of the secondary star. The calcula- 
tions showed that for a fixed orbital period, the braking 
rates dropped by several orders of magnitude as the polar 
field strength increased from ~ 10 to ~ 100 MG, and the 
magnetosphere of the white dwarf engulfed the companion 
star an d trap ped the wind. Subsequently, Li & Wickrama- 
smg he ( J1998D allowed for more general orientations of the 
dipole moment of the secondary star and showed that the 
largest reductions in the braking rate occurred for the case 
where the dipole moment was nearly (within 20°) aligned 
with the spin axis. 



The ROSAT survey flVoges et al. 199$ jVoges et al, 
200C |) re su lted in a large increa se in the number of known 



AM Hers ( [Thomas et al. 1998 ), and allowed a better sta- 
tistical analysis of the period distribution of the magnetic 
and non-magnetic CVs. At long orbital periods, the num- 
bers of AM Hers falls much more rapidly than those of non- 
AM Hers, a difference attributable to a competition between 
accretion torques and magnetic torques which favors mag- 
netic torques at low donor masses, short orbital periods, and 
small orbital separations. As CVs evolve to shorter orbital 
period, then, a growing fraction achieve synchronism and 
evolve from non-AM Her to AM Her systems. Meaningful 
comparisons of AM Her to non-AM Her period distribu- 
tions are therefore limited to short orbital periods, where 
statistics are presumably dominated by the influx of sys- 
tems from longer orbital periods, as opposed to bir th within 
the period interval of interest. Wheatley's ( 1995 ) applica- 
tion of the two-sided Kolmogorov-Smirnov test to systems 
with orbital periods less than 3^2 yielded a 15% probability 
that AM Her and non-magnetic distributions were drawn 
from the same parent population, a di fferen ce he regarded 
as not significant. (Hellier & Naylor (1998), on the other 
hand, found a similar figure for systems with orbital period 
less than 5 h , and regarded it as significant.) We shall see 
below that the differences between AM Her and non-AM 
Her distributions should be most evident within the period 
gap itself, and that a two-sided Kolmogorov-Smirnov test 
applied only to systems within the gap yields a 8.5% proba- 
bility that they belong to the same parent distribution. This 
comparison is still hampered by small numbers, but we re- 
gard the difference as compelling, if not yet conclusive. 

In this paper we reports results of calculations of the or- 
bital evolution of magnetic CVs which allow for a re ductio n 
in magnetic braking following the theory of Li et al. ( 1994a ). 



We use our results to analyse the data on the period distri- 
bution of the magnetic CVs, and more generally to re-assess 
the role played by disrupted magnetic braking in the forma- 
tion of the period gap. 



2 OBSERVATIONAL CONSTRAINTS 

Magnetic field measurements from Zeeman and cyclotron 
spectroscopy provide estimates either of the effective dipolar 
photospheric field, or of magnetic fields at accretion spots 
on the white dwarf surface. We use as our data base for 



AM Hers the list of field m easurements given in (Wick- 
ramasinghe & Ferrario 2000) based on these techniques 



upda ted to include the r ecent field measurement of V884 
Her ( [Schmidt et al. 2001 ). The braking calculations, how- 
ever, require the magnetic moment which is difficult to esti- 
mate from observations. The magnetic moment while being 
proportional to the magnetic field is also proportional to the 
cube of the stellar radius. Unfortunately, the latter, which is 
derived from the mass, is poorly determined. Indeed, even in 
the few eclipsing systems where the orbital inclination is well 
constrained, there is no consensus on the component masses, 
and these are even less well determined in other systems. In 
our calculations, the magnetic moments have been chosen 
to be in the range 10 32 — 10 35 G cm 3 which is consistent 
with the observed effective dipolar field range if allowance is 
made for a spread in white dwarf mass. 

Our list of CVs is drawn from the database of Downes et 



al. (2001) current on 2002 January 25. We have included the 
intermediate polars in our 'non-magnetic' sample because 
to a first approximation, we expect them to evolve under 
the normal braking law. However, it is likely that many IPs 
become synchronised as they evolve toward shorter periods, 
becoming AM Her systems. This possibility needs to be con- 
sidered when comparing observed distributions. 

We show in Figure 1 the distributions of AM Hers and 
other types of CVs as a function of orbital frequency (hr~ ), 
binned into frequency intervals 0.025 hr _1 . In these units the 
canonical 'period gap' corresponds to the frequency range 
0.325 — 0.475 hr _1 . We choose to display the distribution 
in terms of frequency, rather than more conventionally in 
terms of period, in order to spread the short-period peak and 
to make the period gap more apparent. In this figure, and 
throughout the rest of this paper, we will continue to display 
results in terms of orbital frequency for ease of comparison; 
however, we will cast our discussion in the more familiar 
terms of orbital period, and show at the top of each figure 
the corresponding period scale. 

A glance at Figure 1 shows the following characteristics: 

There is an abrupt drop in the number of non-AM Her 
systems at periods above 2.1 hours, which can be considered 
the lower end of the gap. At longer periods, there is a gradual 
rise in their number distribution, with the distribution rising 
rapidly at ~3 hours, which can be considered the upper end 
of the period gap. This plot clearly shows that the gap is 
itself not totally devoid of stars. Rather there is a tendency 
for this canonical period gap to fill in from the higher period 
end. 

The period gap is decidedly less well-defined among the 
AM Her stars. The average number of systems per unit or- 
bital frequency within the 2.1-3 hour period range is smaller 
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Figure 1. Orbital frequency distribution of AM Herculis type 
systems (bottom panel), and all other types of CVs (top panel). 
Filled portions of the histograms correspond to systems with con- 
firmed orbital periods, with the black portions highlighting ob- 
jects in the canonical 2 — 3 h period gap. Open portions of the 
histograms represent orbital periods of additi onal s ystems con- 
sidered uncertain/unreliable by Downes et al. (2001), or derived 
therein from superhump periods. 



than among systems at shorter period (higher frequency), 
but the increase in numbers at longer period (lower fre- 
quency) is much less pronounced than among the non-AM 
Her systems. Moreover, among the magnetic systems, the 
canonical period gap appears to fill from the lower period 
(higher frequency) end, opposite to the behaviour seen for 
the other CVs. The Kolmogorov-Smirnov two-sample statis- 
tic, applied only to those systems of confirmed orbital period 
falling within the gap in Figure 1, yields only a 8.5% prob- 
ability that AM Her systems and non-AM Her systems are 
drawn from the same population. 

We show in Figure 2 the distribution of field strength 
as a function of orbital frequency for the 37 systems for 
which field estimates are available. The fields correspond 
to effective dipolar fields for cases where photospheric Zee- 
man measurements are available, or where cyclotron fields 
are available for both poles. For other systems, the fields 
correspond to the value measured at a single pole. This 
plot shows that there is no significant difference between 
the mean field of systems within the canonical period gap 
and outside. Furthermore, there is no obvious correlation 
between field strength and orbital period for systems within 
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o 




Figure 2. The logarithm of the magnetic field plotted against 
orbital frequency for all AM Hers for which field estimates are 
available. 



3 THE MAGNETIC BRAKING OF 
SYNCHRONISED BINARIES 

3.1 The parametrisation 

The braking rate of synchronised magnetic binaries depends 
on the interaction of the stellar wind with the combined 
magnetospheres of the components stars. It is therefore a 
function of a number of parameters which fall into two basic 
categories: those which specify the magnetic interaction, and 
those which specify the properties of the wind. 

There are four parameters belonging to the first cate- 
gory which enter directly into the calculations of the braking 
rates. These are the masses Mi and M2 of the primary and 
the secondary stars respectively, the radius R2 of the sec- 
ondary, and the magnetic moment fii of the primary star. 
The magnetic moment /12 of the secondary is not a free pa- 
rameter but is calculated on the assumption that prior to 
the star becoming fully convective, the mean surface field 
evolves with angular velocity according to a power law 



Bo = 



UJ 
UJQ 



Br. 







the 2 h — 3 h period region. 



where uj is the angular frequency of rotation. Mestel & 
Spruit (1987) considered the cases p = 1 as being most con- 
sistent with the observations of late type stars, and we ac- 
cordingly adopt p = 1 in our models. However, it should be 
noted that this law does not allow for complexities such as 
dynamo saturation which may play a role in binary evolu- 
tion where rapid rotation rates are encountered as the orbit 
shrinks. 

The magnetic moment of the secondary star 
H2 = O.5B2R2 

is a strong function of its radius i?2- Among the CV donors 
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we wish to model, that radius is itself strongly corre- 
lated with the stellar mass. We therefore parametrise R2 
in terms of a pseudo-radius 7Z2 which removes the lowest- 
orde r p ower-law relationship between R2 and M2 (Warner 
1995^_]riius, 

R2 = fR.2 , 
where 

■JI2 = O.O94R (M 2 /O.O65M ) 13/15 . 

After the star becomes convective we assume that the 
intrinsic mean field is effectively zero. The secondary star 
will have an induced dipole moment of 

fi 2 ^ p\{R2/af 

where a is the orbital separation. The induced dipole mo- 
ment is typically a factor ~ 10 — 100 times lower than the 
intrinsic dipole moment of the primary star, and there will 
be no open flux to produce braking even if a wind persists. 
In this regime we can therefore assume that braking is domi- 
nated by angular momentum loss by gravitational radiation. 

It can be seen from the above, that apart for the wind 
parameters to be discussed in the next section, the braking 
rate can be parametrised terms of Mi, /, /xi, M2. 



3.2 The calculations of the magnetic braking rate 

In the braking of single stars by a magnetised stellar wind, 
both thermal and centrifugal effects can play an important 
role in differ ent mass and rotation regimes. The Mestel & 
Spruit (1987) theory allows a convenient parametrisation of 
these effects. In the context of the CVs, however, due to 
their rapid rotation we can assume that centrifugal effe cts 
dominate over the entire observed range ( |Li et al. 1994a| ; 



et al. 1994b) 



The formulation of the magnetic braking problem for 
centri fugally driven winds has been described in detail by Li 
et al. (1994a). We therefore present only the basic equations 
and our method of solution, and refer the reader to this 
paper for further details. 

The basic assumption underlying our calculations is 
that the dipole moment of the secondary star is oriented 
perpendicular to the orbital plane - the likely orientation 
for a dynamo generated field. The magnetic moment of the 
primary star must then be anti-aligned with that of the sec- 
ondary in the minimum energy configuration (stable equi- 
librium). The magnetic moment of the white dwarf rotates 
through 7r radians as the polarity of the secondary star 
switches during its magnetic cycle. We assume that this ad- 
justment occurs on a time scale that is shorter than the 
magnetic cycle time scale, and that a given system spends 
the bulk of its evolution with the dipoles nearly anti-aligned. 
In this context we note that in the few well studied AM Hers 
there is a pref erence for the dipole axes to be aligned with 
the spin axis (Warner 1995). Changes have been observed 



in the rotational azimuth of the magnetic axis in three of 
AM Her systems on a time scale of ~25 years (see Wickra- 
masinghe & Wu 1991, Wu & Wickramasinghe 1993) but no 
changes in magnetic polarity have so far been reported. 

In what follows we make the explicit assumption that 
the braking rates that we deduce for anti-aligned systems 



are representative of the mean braking rate under which the 
AM Hers evolve. 

We consider a coordinate system in which the secondary 
star is at the origin O, OZ in the direction of the orbital 
angular momentum, and the primary star is in the negative 
y direction at (0,— a,0). The field components in the plane 
OYZ containing the two dipole moments are 



By = 



and 



-3/12Z 



r° rf 



B, 



^2 



+ 



3/cm 



-3^ 



1 + ^hn 

r 5 r 5 



where 

kn = 

where n and r are the radial distances of a field point from 
the primary and the secondary stars respectively. 

The dead zones are assumed to be isothermal with 
sound speed c 2 (d) and to corotate with the secondary star. 
The condition of magnetostatic equilibrium then yields the 
density within this region. For a point with polar coordi- 
nates (r, 8) along a field line starting from the surface of t he 
secondary (r = R, 8 = 80), the density is (Li et al. 1994a) 



P 



po exp 



1 



-Id 



R 



+ - X ld (sin 2 9 



n(R) 



■ P ( — sin 9 — sin i 
* R 



where po at the density at the base of the dead zone. Here 
R 2 lo 2 „ Mia R 2 uj 2 



Xh = 
and 
U = 



p = 



M!+M 2 c 2 (d) 



GM 2 



Ml — 



GMi 



Kd)n(R) 



respectively measure the importance of centrifugal and grav- 
itational effects of the two stars relative to thermal effects 
on the wind. We trace the magnetic field lines by solving the 
coupled differential equations 

dy _ By 
Is ~ ~B 
and 

dz B z 
ds ~ ~B 

where B is the total field. The critical field lines which sepa- 
rate the wind zones from the dead zones are determined by 
the pressure balance condition 

^ = pc 2 M. 

For assumed dea d zone parameters, which we have evaluated 
as in Li et al. ( 1994a| ), these calculations yield the region 
81 < 9 < 92 of open field lines in the plane OYZ which is 
then used to estimate the fractional area $ on the surface 
of open magnetic flux. 

$ depends primarily on the magnetic moment of the 
primary star and of the mass of the secondary star, and only 
weakly on Mi . The dependence of 3> on these two parameters 
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Figure 3. The fraction <E> of open field lines as a function of 
donor star mass M2 and white dwarf magnetic moment m, for 
an adopted white dwarf mass Mi =0.7 Mq and donor radius fac- 
tor / = 1. The shaded region corresponds to complete closure of 
the magnetosphere; in this case, magnetic braking is completely 
suppressed, and systems evolve only through gravitational radia- 
tion. $ is quite insensitive to Mi and /. 



is shown in Figure 3. We note in particular the rapid drop in 
open flux for large fii at fixed M2 which plays a crucial role 
in determining the characteristics of the orbital evolution to 
be discussed in the next section. 



As shown by Li et al. (1994a), the braking rate takes 
a particularly simple form when the wind is centrifugally 
driven. In this limit 



J 



4.63 x lO^dyn cm 

5/3 



Jo 



5 n 

Porb 



4.63 x 10 35 

5 / 3 / B„ \ 10/3 



V0.258/ V0.5/ 



The same theory also provides an expression for Jo which 
involves wind and dead zone parameters. However, the de- 
pendence of these quantities on spectral ty pe is unknown, 
and as discussed by Li & Wickramasing he ( |1998| ), a simple 
scaling from solar values is inappropriate. Following previ- 
ous investigators, we have accordingly chosen to normalise 
the braking rate so as to reproduce the observed period gap 
in the limit of zero magnetic field. This normalisation corre- 
sponds to a braking rate of 4.63 x 10 35 dyn cm for a Roche 
lobe filling secondary star in a non-magnetic (n\ = 0) CV 
at an orbital period of 5 h (Rappaport et al. 1983). 

The mass of the white dwarf and its magnetic moment 
are essentially free parameters in our calculations as are the 
mass of the secondary and the factor / which determines 
the radius of the donor star. We have constructed a grid of 
magnetic braking rates as a function of Mi , /, /ii , M2 , and 
used these to calculate orbital evolution for various initial 
component masses. 
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Figure 4. Mass transfer rates as functions of orbital frequency 
for different degrees of systemic mass loss. The parameter f3{= 
1.0,0.5,0.0,-0.5,-1.1) is the fraction of mass lost by the donor 
which is retained by the white dwarf. Where (3 = 1, mass transfer 
is completely conservative, while for (3 = — 1 the white dwarf is 
presumed to expel twice the mass it accretes over a nova outburst 
cycle. Our assumed value (/3 = 0) is emphasized. The initial donor 
mass is fixed at M^i = 0.6 Mq, and the initial white dwarf mass 
adjusted to yield the same total mass (Mi + M2 = 1.268 Mq) at 
the upper edge of the period gap in all of the sequences shown. 



4 EVOLUTIONARY CALCULATIONS 

We have explored the evolutionary consequences of atten- 
uated magnetic braking usi ng th e bipolytrope code devel- 
oped by Rappaport et a l. ( 198J ) , as modified by Nelson, 
Rappaport & Joss ( 1993 ). This code idealizes the convec- 
tive envelope-radiative core structure of near-main sequence 
stars of arbitrary degeneracy as n = 3/2 polytropic en- 
velopes fitted to n = 3 polytropic cores, and reproduces the 
most important features of mass-losing lower main sequence 
stars with considerable fidelity. Except as noted below, we 
assume solar composition for the donor star throu ghout , 
with photospheric opacities adopted from Alex ander (197E; 



198C| ) after Rappaport, Joss & Webbink ( |1982j ). 

For purposes of computing mass transfer rates, we fit 
isothermal atmospheres to the bipolytropes, and write the 
mass transfer rate explicitly in terms of the difference be- 
tween the fitting point radius for this atmosphere and the 
Roche lobe radius of the donor star. However, in the calcula- 
tions presented here, that difference never exceeds 2 x 10~ 3 
of the stellar radius, so we could have achieved the same re- 
sults for all practical purposes by limiting the stellar radius 
to the Roche lobe radius as a surface boundary condition. 

We assume that the white dwarf mass remains constant 
throughout the evolution of our model binaries, with all 
mass lost by the hydrogen-rich donor lost from the system 
in nova explosions. The ejected mass is assumed to carry 
off a specific angular momentum equal to the orbital angu- 
lar momentum per unit mass of the white dwarf. Observa- 
tionally, nova outbursts are found in non-AM Her systems 
and AM Her systems alike {viz. V1500 Cyg). The heavy- 
element enhaiu^ernen^t^ypical of the ejecta of well-studied 
novae (Gehrz et al. 1998) are unmistakable evidence that 
the underlying white dwarfs do not generally grow in mass 
during the course of evolution, but may in fact suffer net 
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Figure 5. Mass transfer rates as functions of orbital frequency 
for systems differing only in initial donor star mass (M%i = 
1.0,0.8,0.6,0.4 Mq). The emphasized curve corresponds to the 
same initial conditions (Mu = 1.0 Mq, M^i = 0.6 Mq, /ii=0, 
f3 = 0, and solar mctallicity, [Fe/H]=0, throughout Figures 4-8. 
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Figure 6. Mass transfer rates as functions of orbital frequency 
for systems differing only in initial white dwarf mass (Mu = 
1.4,1.2,1.0,0.8,0.6,0.4 Mq). The emphasized curve is common 
to Figures 4-8. 



erosion through nova outbursts. Provided that this erosion 
is not so severe as to drive the binary toward dynamical in- 
stability (roughly, that dMi/dM 2 < 3/2), that the binary 
in question has not arrived at the onset its cataclysmic vari- 
able phase just above the period gap, and that the angular 
momentum content of nova ejecta has not been seriously un- 
derestimated, then the width and position of the period gap 
are insensitive to our assumption of constant white dwarf 
mass. Figure 4 illustrates this point, by comparing a fam- 
ily of evolutionary tracks for nonmagnetic CVs, differing in 
the extent of systemic mass loss, but converging to identical 
total masses at the upper edge of the period gap. 

Figure 5 underscores the point implicit in Figure 4: for 
a given magnetic braking rate, the position and width of the 
period gap depends principally on the component masses of 
the binary as it ceases magnetic braking, and not on its evo- 
lutionary history in reaching that point. Since the disrupted 
magnetic braking model for the period gap associates the 
cessation of magnetic braking with the donor star becoming 
fully convective, the principal factors governing the position 
and width of the period gap (again, for a given magnetic 
braking rate) are (i) the mass of the white dwarf (or, equiv- 
alently, the mass ratio), and (ii) the chemical composition 
of the donor star. 

The mass of the white dwarf influences the period gap 
through its role in orbital dynamics. Smaller white dwarf 
masses demand larger mass transfer rates to moderate the 
rate of Roche lobe contraction induced by angular momen- 
tum losses to magnetic braking. These higher mass trans- 
fer rates in turn drive the donor star further from thermal 
equilibrium, expanding the period gap, as seen in Figure 6. 
However, the effect is only strong for binary mass ratios near 
the limit of dynamical stability for fully convective donors 
((Mi/M2) C rit = 3/2), and cannot by itself account for CVs 
found in the gap. 

Figure 7 illustrates the effect of varying metallicity of 
the donor star on the position and width of the period 
gap. Metal-deficient donors are systematically bluer and 
more luminous than their solar-composition counterparts 
of the same mass. Their higher luminosities endow them 
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Figure 7. Mass transfer rates as functions of orbital frequency for 

systems differing only in donor star metallicity ([Fe/H]= 0.0, 

0.5, 1.0, 1.5, 2.0). The emphasized curve is common 

to Figures 4-8. 

with shorter thermal time scales, and they therefore adhere 
more closely to the thermal equilibrium main sequence be- 
fore turning fully convective than do solar-composition stars. 
As seen in Figure 6, then, metal-deficient stars encroach into 
the period gap from the long-period (low-frequency) end. We 
tentatively suggest, therefore, that the small number of non- 
AM Her-type CVs spilling into the period gap at long period 
may be relatively metal-poor. Metal-deficient main sequence 
stars are also physically smaller at low masses than solar 
composition stars, and so the short-period (high-frequency) 
boundary of the period gap migrates toward shorter periods 
(higher frequencies). 

The preceding considerations demonstrate that, within 
the context of the disrupted magnetic braking model, the 
period gap is a robust entity. One can account for small num- 
bers of non-AM Her CVs in the gap by appealing to newly- 
formed systems in which the donor stars are already fully 
convective main sequence stars with masses in the 0.2-0.3 
Mq range, or by appealing to infiltration by metal-deficient 
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Figure 8. Mass transfer rates as functions of orbital frequency 
for systems differing only in white dwarf magnetic moment (fix = 
0, 10 33 , 10 34 , 10 35 G cm 3 ). The emphasized curve is common to 
Figures 4-8. 



donors at the long-period end of the gap. However, in light 
of the robustness of the period gap, the difference in orbital 
period distributions between AM Her systems and non-AM 
Her systems seen in Figure 1 strongly suggests that the AM 
Her systems deviate from the disrupted braking law which 
seems to apply to all other CVs. 

As described in Section 3, in synchronous rotators (AM 
Her binaries), the dipole field of the white dwarf can close 
field lines from the donor star, suppressing the stellar wind 
on otherwise open field lines, and thereby suppress angu- 
lar momentum losses. We have incorporated the attenuated 
magnetic braking model described there into the evolution- 
ary model described above, under the assumption that the 
white dwarf magnetic dipole moment, like its mass, remains 
constant in time. Our model sequences assume that the 
white dwarf is locked in synchronism with the binary orbit 
from the very onset of mass transfer, i.e., that the system is 
of AM Her type throughout its evolution. 

We show in Figure 8 the effect of increasing magnetic 
moment, /lii, of the white dwarf primary on the evolution 
of a typical AM Her system with initial component masses 
Mi = 1.0 Mq for the white dwarf and M 2 = 0.6 Mq for 
the donor secondary. As the braking torque is suppressed 
with increasing fii , mass transfer rates above the period gap 
subside, and the donor star is no longer driven so far from 
thermal equilibrium before it becomes fully convective at the 
top (lower frequency edge) of the period gap. As the donor 
star is no longer so severely inflated in comparison with its 
thermal-equilibrium main sequence counterpart, so its con- 
traction following cessation of the magnetic braking driv- 
ing mass transfer is also less severe. Mass transfer resumes 
at longer orbital period (lower orbital frequency) than in 
the disrupted magnetic braking model for non-AM Her-type 
CVs because the donor is no longer so severely undermas- 
sive for its orbital period, and the period gap fills from its 
shorter period (higher frequency) end. As seen in Figure 5, 
this filling of the period gap is progressively greater as 
increases, and magnetic braking prior to the secondary's loss 
of its dipole field is suppressed. For sufficiently high values 
of fii (sec \x\ = 10 35 G cm 3 in Figure 5), magnetic brak- 
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Figure 9. Upper and lower limits to the period gap, as functions 
of white dwarf mass M\ (as labeled, in units of Mq) and white 
dwarf magnetic moment The gap boundaries depend on 
only for initial masses M^i < 0.4 Mq. fi = is assumed. 



ing becomes inconsequential, and the system evolves more 
or less continuously through the period gap (the momentary 
dip in mass transfer rate at P « 3 h owing to an inflection 
point in the thermal equilibrium mass-radius relation at the 
point stars become fully convective). 

Figure 9 illustrates the dependence of the upper and 
lower edges of the period gap on white dwarf mass Mi and 
magnetic moment . We define the upper (lower-frequency) 
edge of the period gap as the point at which the donor 
star becomes fully convective. In those cases where mag- 
netic braking remains strong enough to drive the donor out 
of thermal equilibrium, this point indeed marks a local min- 
imum in the orbital period, which increases momentarily (if 
only by a very small increment) as the mass transfer rate 
subsides. We take the lower (higher-frequency) edge of the 
period gap as the period maximum (frequency minimum) 
during the thermal time scale mass transfer episode visi- 
ble as a slim vertical loop in Figure 5. For sufficiently large 
white dwarf magnetic moments, where magnetic braking is 
nearly completely suppressed, this loop disappears, and we 
then define the lower edge of the period gap as the local 
maximum in M2. 

It is evident from Figure 9 that, for white dwarf mo- 
ments fii 10 33 G cm 3 , the gap edges asymptotically ap- 
proach those of the non-AM Her systems. At the other ex- 
treme, for white dwarf moments fii ^ 4 x 10 34 G cm 3 , the 
gap is reduced to a momentary hiccup as the donor becomes 
fully convective. For intermediate white dwarf moments, the 



© 2001 RAS, MNRAS 000, ??-?? 



R.F. Webbink & D.T. Wickramasinghe 



3.0 



m 



o 



8.5 



8.0 



7.5 



7.0 



6.5 



1.2 M„ 



RX J050 1-0359 




(hr) 



2.0 



0.35 



0.40 



0.45 



0.50 



(hr"') 



Figure 10. Observed field strengths of AM Her stars as functions 
of orbital frequency, in the region near the period gap. The con- 
tinuous curves correspond to model calculations (cf. Fig. 9) with 
differing white dwarf mass Mi (= 1,2, 1.0,0.8,0.6 Mq), and en- 
close the period gap corresponding to the respective value of M\ . 



period gap fills from its short-period (high-frequency) end, 
as indeed observations suggest (see Figure 1). A clear pre- 
diction of our model is that, as more AM Her systems are 
discovered within the period gap, a clear trend of increasing 
magnetic moment with increasing orbital period (decreasing 
orbital frequency) should emerge. 

To facilitate comparison with observed AM Her stars in 
the period gap, we have used the Nauenberg (1972) mass- 
radius relation for cold white dwarfs to convert the dipole 
moments of Figure 9 to polar fields. The region immedi- 
ately surrounding the period gap is shown in Figure 10. Un- 
fortunately, empirical masses of white dwarfs in AM Her 
binaries are generally very poorly-constrained. However, to 
the extent that our model results can be taken at face 
value, the positions of observed systems in Figure 10 can 
be used to set upper limits to their white dwarf masses, 
since for given dipole moment, the polar field strength de- 
creases rapidly with increasing mass. In this figure, AM Her 
binaries must fall outside the gap corresponding to their 
white dwarf masses, thus implying that QS Tel, V381 Vel, 
V516 Pup, and UW Pic must contain white dwarfs of mass 
Mi ^ 1.1 Mq. The reader is cautioned, however, that this 
constraint is quite sensitive to the precise location of the 
short-period (high-frequency) edge of the canonical period 
gap. 



5 DISCUSSION AND CONCLUSIONS 

The AM Her stars, unlike cataclysmic variables generally, 
have mostly been discovered through all-sky surveys, where 
they appear as strong X-ray emitters. We can therefore ap- 



°o 10 5 4 



1.5 



1.0 



Spherical 
source 
distribution 




Figure 11. Relative discovery probabilities for AM Her stars as 
a function of orbital frequency. Only results for our fiducial initial 
model (Mi = 1.0 Mq, M 2 i = 0.6 Mq, = 0, solar metallicity) 
are shown, for different values of white dwarf magnetic moment 
fil, as in Figure 8. A constant formation rate over the age of 
the galactic disc (10 10 yr) is assumed. The nonmagnetic case is 
hatched, with the hatching sloping upward in the direction of 
evolution. The vertical dashed segments mark singularities which 
arise in the distributions shown where !> or t, changes sign. Such 
singularities occur in pairs at the upper and lower edges of the 
period gap, but only those at the lower (higher-frequency) edge 
of the gap are resolvable here. 



proximate the selection effects governing their discovery by 
treating them as an X-ray-brightness limited sample, with 
X-ray luminosities proportional to the mass transfer rate 
\M2\- To the extent that X-ray extinction is negligible, the 
distance to which systems are sampled is then proportional 
to IM2I 1 - If the threshold distance is small compared with 
the scale height of AM Her systems with respect to the 
galactic plane, then the apparent distribution should be 
roughly isotropic, and the number of systems sampled varies 
as IM2I 3 , should the distance be large compared with the 
scale height (but still small compared with the radial scale 
length of the galactic disk) the number of systems varies as 
IM2I itself. The number of systems in any orbital frequency 
interval is then simply proportional to the lifetime of the 
system in that interval, i.e., to v~^. As a point of reference, 
the mean value of sin \b\ of the AM Her systems in Figure 
1 is 0.51, indistinguishable from that of a spherical source 
distribution ((sin |6|) = 0.50). 

Figure 11 illustrates the relative discovery probabilities 
arising from this simple model of selection effects, in the 
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two limiting cases regarding their spatial distribution, for 
our fiducial evolutionary sequence. In reality, we expect CVs 
to form with a range of component masses, and this will 
affect the overall form of the true distribution. In particular, 
systems with a given initial secondary mass, M21, contribute 
to the distribution only at orbital frequencies ^ or b above 
the frequency at which the CV first reaches interaction (cf. 
Figure 5). Likewise, magnetic CVs will contribute to the 
AM Her distribution only at orbital frequencies above those 
at which they first reach synchronism. Nevertheless, as seen 
above, the properties of individual systems at the upper and 
lower edges of the period gap, and the positions of those 
edges, are insensitive to their initial conditions. 

One feature immediately evident from Figure 11 is that 
the number density per unit orbital frequency of CVs at 
the upper (low-frequency) edge of the period gap is decid- 
edly greater than at its lower (high-frequency) edge. The 
ratio of these two number densities, as defined at the very 
edges of the period gap, amounts to a factor of 5.4 for non- 
AM Her systems with a spherical source distribution (1.9 
in a planar distribution). Among the AM Her systems of 
increasing magnetic moment fix, the amount by which this 
ratio exceeds unity is roughly proportional to the width of 
the period gap appropriate to that value of fix- 

Reference to Figure 1 shows that the non-AM Her sys- 
tems have a much higher number density of systems at the 
upper edge of the period gap than at its lower edge, in at 
least qualitative agreement with Figure 11. However, among 
AM Her systems themselves, this ratio of number densities 
is plainly less than unity. The period gap is for them only 
partially filled, and the marked increase in number density 
beyond the lower (high-frequency) edge of the gap should 
be matched by an equally marked increase at the upper 
(low-frequency) edge. The clear implication is that more AM 
Her-type systems emerge from the period gap than enter it. 
Thus, a significant fraction of AM Her systems below the 
period gap must have entered the gap in a non-AM Her 
(asynchronous) state, i.e., as intermediate polars. As asyn- 
chronous systems, they suffered full magnetic braking above 
the gap, and were driven out of thermal equilibrium to the 
same extent as other non-AM Her systems. They therefore 
remain detached until resuming mass transfer at the lower 
edge of the gap as defined by other non-AM Her systems, in 
transit becoming synchronously locked, but failing to appear 
either as CVs or AM Her systems within the gap. Systems 
which achieve synchronism while crossing the period gap 
thus contribute to the AM Her population only at periods 
below the gap. 

Comparison of Figure 11 with Figure 1 also shows that 
the theoretical period minimum of AM Her stars lies at sig- 
nificantly shorter period (higher frequency) than observed. 
This is a problem pervading CV evolution generally, since 
both AM Her systems and non-AM Her systems are driven 
through this per iod m inimum by gravitational radiation. 
Kolb & Baraffe (1999) discuss a number of possible expla- 



netic field should suppress a magnetic stellar wind entirely 
among AM Her systems. The fact that both AM Her bi- 
naries and non-AM Her binaries are subject to the same 
minimum period (Figure 1) then suggests that stellar winds 
are not responsible for this discrepancy. 

We summarise our principal conclusions as follows: 

(a) The width of the period gap, and the position of 
its lower boundary, is a function of the magnetic moment 
of the primary star. For sufficiently high magnetic moments 
(fix il 4x 10 34 G cm 3 ), the gap disappears, and mass transfer 
is driven primarily by gravitational radiation. 

(b) Observations of AM Her-type systems show that the 
period gap is less pronounced than that of CVs belonging to 
all other classes. In particular, the AM Hers appear to fill 
the gap from the short-period (~2 h ) end, behaviour which 
supports the reduced braking theory outlined in section 3. 
In contrast, non-AM Her systems encroach on the period 
gap from its long-period edge, behaviour which we suggest 
may reflect lower metallicities of these particular systems. 

(c) The differences between the distributions of AM Her 
binaries and non-AM Her binaries within the period gap 
in turn provide strong prima facie evidence that magnetic 
stellar winds are indeed the dominant angular momentum 
loss mechanism among non-AM Her binaries. 

(d) The abundance of AM Hers immediately below the 
period gap, in comparison with their numbers immediately 
above the period, requires that a significant proportion of 
the AM Hers found below ~2 h entered the period gap as 
IPs, and become synchronised only as mass transfer ceased 
within the gap. The smaller accretion torques, and increased 
magnetic interaction (due to the reduction in orbital separa- 
tion) allowed these systems to remain synchronised as they 
emerge from the gap. 
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